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Abstract 

We observe a condensate forming in the interior of a black hole (BH) during numerical 
simulations of gravitational collapse of a massless charged (complex) scalar field. The 
magnitude of the scalar field in the interior tends to a non-zero constant; spontaneous 
breaking of gauge symmetry occurs and a condensate forms. This phenomena occurs in the 
presence of a BH without the standard symmetry breaking quartic potential; the breaking 
occurs via the dynamics of the system itself. We also observe that the scalar field in the 
interior rotates in the complex plane and show that it matches numerically the electric 
potential to within 1%. That a charged scalar condensate can form near the horizon of a 
black hole in the Abelian Higgs model without the standard symmetry breaking potential 
had previously been shown analytically in an explicit model involving a massive scalar field 
in an AdSi background. Our numerical simulation lends strong support to this finding, 
although in our case the scalar field is massless and the spacetime is asymptotically flat. 


1 Introduction 

In this work, we consider the gravitational collapse of a complex (charged) massless scalar 
field to a charged black hole with asymptotically flat spacetime. Our model corresponds to 
scalar electrodynamics coupled to gravity (with the gauge field coupled to the scalar field in 
the usual fashion via the covariant derivative.) However, the standard symmetry breaking 
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potential V(cj)) = n 2 q i 2 + A< f) 4 (with fi 2 < 0 and A > 0) is absent. Despite this absence, we 
observe numerically that the scalar field in the interior tends to a non-zero constant in contrast 
to the collapse of uncharged matter. Spontaneous breaking of the gauge symmetry occurs and 
a condensate forms. The breaking occurs via the dynamics of the system itself. In other words, 
in the presence of the charged black hole, it is energetically favorable for the magnitude of the 
scalar field to be non-zero in the interior. We also observe that the Ricci scalar at late times 
has a non-zero constant value inside confirming that a condensate of constant energy density 
has indeed formed in the interior. 

The scalar held is actually observed in the interior to rotate with time in the complex plane. 
From theoretical considerations, the angular speed uj should be equal to the “electric potential” 
a in the interior. We show that the numerical value for u matches a to within 1%. Previous 
analytical work had already noted that despite the absence of the symmetry breaking potential 
a charged scalar condensate can form near the horizon of a BH [i][2] . This was shown 
explicitly in a model with negative cosmological contant and massive scalar held in an AdS^ 
background [T]. In our case, the spacetime is asymptotically hat and the scalar held is massless. 
Despite this difference, we observe numerically that a charged scalar condensate forms as in (I]. 
In hat spacetime, it was already shown that spontaneous symmetry breaking (SSB) can occur 
without including explicitly the standard symmetry breaking potential [3j. The authors discuss 
a system at finite charge density (i.e. finite chemical potential) where SSB occurs but is absent 
when the charge density is zero. They find that states of definite charge that are local minima 
of the energy correspond to the complex scalar held rotating with time with the angular speed 
uj identified with the “chemical potential”. As already stated above, the rotation of the scalar 
held is precisely what we also observe, with the angular speed uj identihed in our case with 
the “electric potential”. One of the boundary conditions we impose on the charged collapse is 
that the electric potential is held hxed asymptotically. Thermodynamically, this corresponds 
to the grand canonical ensemble and the electric potential is then identihed as the chemical 
potential. This clearly shows the connection between our black hole case and the hat spacetime 
case discussed in j3]. 

The code for the numerical simulation, the initial conditions and the equations of motion 
governing the charged collapse are the same as those found in our previous work [1] , where the 
focus was on BH thermodynamics. Therefore, we will only provide enough information and 
equations in the next section so that the reader can follow in later sections the presentation of 
results pertaining to the formation of the condensate. For the upcoming section, the reader can 
always refer to j4] if they would like to see derivations or more details. We work in geometric 
units where G = c= 1. We also set Coulomb’s constant to unity. In these units, radius, mass, 
energy, time and charge have units of length. The coupling constant e has units of inverse 
length. It’s value is arbitrary as there is no experimental value for it in our model. We set it 
equal to unity and leave the scale unspecified (one can work in centimetres, metres, etc.). One 
unit of time is the amount of time it takes light to travel a radial distance of one unit in flat 
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spacetimc0. The metric has Lorentzian signature (—, +, +, +). Spacetime indices are in Greek 
and run from 0 to 3. 


2 Evolution and constraint equations 

2.1 Metric in isotropic coordinates 


During the numerical simulation, the metric is time-dependent and spherically symmetric. In 
isotropic coordinates it takes the form 

ds 2 = -N(r,t) 2 dt 2 + iP(r,t) 4 {dr 2 + r 2 dtt 2 ). (1) 


The function t ) is often referred to as the conformal factor and the function N(r, t ) is called 
the lapse function. We assume asymptotic flatness so that N and are both unity at infinity. 
Note that the areal radius is given by ) 2 r. The metric © is convenient to use in numerical 
simulations of spherically symmetric collapse because the metric function is finite at the 
horizon [Mo]. It is also useful for the study of BH thermodynamics because the coordinate 
time t coincides with the time measured by a clock at rest at infinity; the Lagrangian can then 
be identified with the negative of the free energy mm- The Reissnner-Nordstrom (RN) BH 
in isotropic coordinates is given by mum 


ds 2 = 


1 - 


M 2 -Q 2 \2 
4r 2 


M M 2 -Q 2 \i 
1 + — ^ 4^2 


dt 2 + 


M , M 2 -Q 2 
1 + — ^ 4^2 


i 2 


[dr 2 + r 2 dn 2 ^j . (2) 


The outer horizon in isotropic coordinates is situated at r + = \f M 2 — Q 2 /2. At this location, 
the lapse function is zero (lV(r+) = 0) and the conformal factor is finite and given by i />(r+) = 
[2 + 2 M/{M 2 — Q 2 ) 1 / 2 ] 1 / 2 . For the Schwarzschild case (Q = 0), V’( r +) = 2 in agreement 
with the uncharged case [9]. The metric ((2]) is valid only in the exterior static region. The 
coordinate r in ([2]) is valid in the range r > r + (the region r < r+ does not correspond to the 
interior but a second covering of the exterior region). The metric for the interior of the RN 
BH is nonstationary |101ll3j . During our numerical simulation, the metric functions N(r,t ) 
and il>{r,t) appearing in ([T|) are nonstationary in the region r<r + , reflecting the true nature 
of the interior spacetime. It is only in the exterior region r > r+ that the metric © approaches 
the static form ([2]) at late times of the collapse. As pointed out in [3], we do not observe the 

x For example, if the scale were specified to be 3 x 10 5 km, then one unit of time would correspond to one 
second. It can be readily checked that the equations of motion do not depend on the choice of scale for e i.e. 
they are invariant under r — > br, t — > bt and e —> e/6 where 6 is a positive constant. 
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second (inner) horizon or timelike singularity associated with the RN BH. We observe only 
one horizon at r+ (where N(r + ) = 0) and a spacelike singularity as in the Schwarzschild case. 
This is in accord with the findings of previous numerical work on charged collapse H5H32I. The 
inner horizon of the RN BH is an artifact of exact staticity (and exact spherical symmetry) [HI 
and it has been known since pioneering work in the 90’s [18] . that the inner (Cauchy) horizon 
is unstable to perturbations. 


2.2 Equations of motion for matter 


For matter, we consider a complex (charged) scalar field coupled to an electromagnetic field 
A/j,. The matter Lagrangian density C m has a local f7(l) gauge symmetry and is given by 

An = ~ (X;m + ieA »x)g fiV {x-u ~ ieA„x) ~ ( 3 ) 

where a semi-colon denotes covariant differentiation evaluated with metric a bar denotes 
complex conjugation and = A V]IM — A^-^ is the electromagnetic field tensor. Spherical 
symmetry and gauge freedom reduces the number of gauge components from four to one: only 
At = Aq is non-zero. For simplicity we label it a i.e. a = At. The matter fields are therefore 
X = x( r T) an d a = a(r,t). Lagrange’s equations of motion for matter are 


V, 


dC r 

l dc l( 


d£rr. 

dq 


= 0 


where q is a generic field. 

Equations of motion for scalar field, x 
We define the quantity 


ip® 

p ^ iv (*■ + ieax ' 


The evolution equations for p and x respectively are [3] 



d r (Nfi 2 r 2 x') — ieap. 


and 


N 



ieaip 6 x\ 

) ' 


(4) 

(5) 

( 6 ) 

(7) 


Equations of motion for gauge field a 
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The constraint equation for the gauge field a (one can view it as Poisson’s equation) is given 
by 



A p 2 a'r 2 \^ 


2me (XP-XP) ■ 


( 8 ) 


Defining the function g as 


its evolution is given by 


9 = 


a'ip 2 

~JT 


g = 2meNip 2 (xx ~ Xx') ■ 


(9) 

( 10 ) 


The equation governing the gauge field a is not an evolution equation but an ordinary differ¬ 
ential equation: a' = g N/ip 2 . 


2.3 Gravitational sector 


Stress-energy tensor 

The stress-energy tensor T^ v appearing in Einstein’s field equations can be calculated from the 
matter Lagrangian ([3]). Its non-zero components are 


Tee = 
T rr = 
T tt = 


1 (PP 


2\^~ XX + 


1 


- I PP . 
n Ts + X X ~ 


Ip 8 


— I^ + 

- 1 -12 T 


2 

N 2 

~2 


x'x 1 


ip 


ip 4 



—) 

Amp 4 J 



T h = Tee sin 2 6 
Trt = ^{xP + XP) 


( 11 ) 


Field equations 

Einstein’s field equations are given by G^ = k 2 T^ u where k 2 = 8ir G = 8 tt and G^ is the 
Einstein tensor evaluated with metric (JT]). The G rr equation yields evolution equations for 
the conformal factor ip while Gtt and G r t yield constraint equations. Gee yields an ordinary 
differential equation for the lapse function N (G^ yields the same equation). 


Define 



( 12 ) 
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K is the negative of the trace of the extrinsic curvature for spacelike hypersurfaces at constant 
time t. Its evolution is given by 


k _ K 2 Qip' Pip' 1\ 3N' (2ip' 1\ 
N 2 ip 5 \ ip + r) Nijfi \ ip r) 

(Y / xY pp g 2 \ 

4 y 1p A ip 12 47T'0 8 y 


The energy constraint is given by 


4 

ip^r 


(2 ip 1 + rY) 


2 


xY pp g 2 \ 

1p A ^,12 4ti-^) 8 y 


while the momentum constraint is given by 


AT 2 

T 


(13) 


(14) 


A'' 

T 


4^6 (x'p + x'p) • 


The ordinary differential equation governing the lapse function N is 

2r ^ ( W_\ « 2 xY , «V 

ip 2 r l K rip 3 J N r \rip A ) ip A 4:iri/j 8 


(15) 


(16) 


The evolution equation (1131) for K can be made more numerically stable by combining it with 
(fl4l) to obtain 


k _ k 2 1\ 3 N' 

N ip 5 \ ip r) N'ip A 
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rip 5 


(2ip' + rip"^ 


3k 2 g 2 
Snip 8 



(17) 


To summarize, the evolution equations for the matter sector are ([7j) and (|5|) for the pair (x,p) 
and m for g respectively. The evolution equations for the gravitational sector are (1121) and 
(fl7l) for the pair {ip, K). The gauge held a and lapse function N obey the ordinary differential 
equations ([9]) and (fT6l) respectively. There is one constraint equation in the matter sector, 
namely “Poisson’s” equation (j8|) . In the gravitational sector, the energy and momentum con¬ 
straint are given by equations (1141) and (1151) respectively. Once the initial states and boundary 
conditions are specified, the evolution of the fields is unique. 
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2.4 Initial state and black hole formation 


The black hole forms when the lapse function N (r, t ) crosses zero for the first time (this occurs 
at approximately t=8.4 in our simulation). The horizon grows with time i.e. the coordinate 
radius r q where N = 0 grows with time and approaches at late times a value of ro = 0.344 
(see Fig. [I]). Initially, the lapse function is everywhere positive (see Fig. [2]) so that there is no 
black hole at the start of the collapse. 

Let x\ and X 2 be the real and imaginary part of the complex scalar field x- We begin by 
choosing the initial configuration for %i, X2 > Xi and X 2 - We choose a kink shape given by 
Xi = X 2 = A2 (tanh (Ai — r) + 1). Their time-derivatives are chosen to be xi = ~X 2 = 

(this ensures that the momentum constraint is trivially satisfied ). One is free to choose the 
values of the different A,, and these in turn determine the values of the conserved mass M and 
charge Q. For this run, we chose values of Ai = 2, A2 = 0.11 and A3 = 0.032. The initial states 
for the metric functions and N and for the gauge field a are then obtained by solving three 
coupled second order differential equations: the energy constraint (1141) . the ordinary differential 
equation (1161) and “Poisson’s” equation ([8]). The initial state for N, ij} and a are plotted in 
Fig. EJ 


3 Numerical results: condensate formation and rotation of 
scalar field 

3.1 Formation of a condensate 

The magnitude of the complex (charged) scalar field at different times of the collapse is plotted 
in Fig. [3] for an initial kink shape. Initially, the field decreases in the interior but later increases 
to a constant value inside. At late times in the collapse process, the magnitude tends towards 
zero in the exterior region but has a non-zero constant value inside, signaling breaking of the 
Abelian gauge symmetry as was originally found analytically in . A condensate forms in 
the interior. This is in sharp contrast to the uncharged case where the value of the (real) 
scalar field tends to zero in the interior at late times [9]. The Ricci scalar is plotted in Fig. EJ 
Note that at late times it has a non-zero constant value inside. By Einstein’s field equations, 
the Ricci scalar is proportional to the negative of the trace of the energy momentum tensor. 
This implies there is a non-zero constant energy density in the interior acting much like a 
cosmological constant. The non-zero constant value of the Ricci scalar in the interior confirms 
that a condensate of constant energy density has indeed formed in the interioid. 

As pointed out in [4], at late times, all aspects of the exterior RN BH are reproduced in the 

2 There is also a Ricci singularity in the interior which is discussed in detail in [8]. 


7 




numerical simulation. In particular, in the exterior region, the metric functions if> and N are 
basically static at late times and match almost exactly the RN BH metric (|2|). The scalar field 
tends towards zero in the exterior region; there is “no scalar-hair” in agreeement with the “no 
scalar-hair” conjecture for spherically symmetric, asymptotically flat static black holes m- 
The solution in pLJ has “scalar-hair” but the spacetime in that case is not asymptotically flat 
but asymptotically AdS 4 and hence does not violate the “no scalar-hair” conjecture. 

3.2 Rotation of the scalar field 

In the interior, we observe that the scalar field rotates with time in the complex plane (see 
Fig. [5]). This is in accord with the fact that D^x = + approaches zero at late times 

in the interior (this condition minimizes the “kinetic term” for the scalar field). This yields 
two equations: dtx + zeay = 0 and d r x = 0 (recall that the only non-zero component of the 
gauge field A M is A t = a). The gauge field a (which can be viewed as the “electric potential”) is 
set to zero at asymptotic infinity (same gauge choice as in |l|jE At late times, it approaches a 
non-zero constant at the horizon and the interior (see Fig. [6]). This implies that the scalar field 
in the interior is given by x = Re~ ieat = Re~ iat where R is the magnitude of the scalar field 
(e = 1 is the value used in our numerical simulation). The scalar field therefore rotates with 
angular speed equal to the potential a at the horizon. The coordinate t in metric (JTj) remains 
timelike in the interior; it does not switch to being spacelike as in standard coordinates and 
the metric functions t/i(r, t) and N(r,t ) are continuous and finite across the horizon. Therefore, 
the scalar field is truly rotating with time with angular speed a in the interior. The numerical 
data corroborates this. In Fig. [7]we plot the rotation speed uj of the scalar field and the average 
value of a in the interior starting at early times. Note that they converge at late times. We 
also made a table of the data between t = 21 and t = 23.5 (see table [[]). The value of u 
and the average value a of the electric potential agree to within less than 1%. As noted in 
the introduction, SSB with complex scalar field rotating with time had also been found in flat 
spacetime in the absence of an explicit standard symmetry breaking potential [3]. 

3.3 Energy conditions 

It is common for scalar fields to violate some energy conditions. Here we will check the 
weak energy condition (WEC) and null energy condition (NEC) and we will see that they 
are both obeyed during the collapse process. The WEC requires that t^t u > 0 for any 
timelike vector t^. We choose the timelike vector = (1,0,0,0) so that the WEC for that 
vector becomes T tt > 0. We plot T a {r,t) at different times in Fig. 0 It is nonnegative at 
all times and the WEC is obeyed. The NEC requires that T^ u > 0 for any null vector 


3 This boundary condition corresponds thermodynamically to the grand canonical ensemble [3j. 



Time 

20.999-21 

21.499-21.5 

21.999-22 

22.499-22.5 

22.999-23 

23.499-23.5 

OJ 

a 

% diff. 

0.18570 

0.18726 

0.83307 

0.18287 

0.18428 

0.76514 

0.18022 

0.18149 

0.69976 

0.17776 

0.17889 

0.63167 

0.17549 

0.17649 

0.56660 

0.17340 

0.17429 

0.510643 


Table 1: The angular speed uj at which the complex scalar field rotates is compared to the 
average value a of the electric potential at a point in the interior (at r = 0.25). The difference 
is less than 1 %. 

P 1 . We choose the null vector £ IA = (1 /N, 1 /,0,0) so that the NEC for that vector becomes 
P = Tu/N 2 + T rr /ili 4 + 2T rt /(N i/j 2 ) > 0. We plot P(r,t) at different times in figures l9l and flOl 
(Fig. [10] is a magnified version of Fig. [9]). One can see in Fig. [10] that P is nonnegative at all 
times and the NEC is obeyed. We therefore have that the NEC and WEC are both obeyed. 


4 Conclusion 

We have shown that numerical simulations of charged collapse reveal that a condensate can 
form in the interior of a charged BH in scalar electrodynamics. Breaking of the Abelian 
gauge symmetry occurs in the presence of the charged BH in the absence of the standard 
symmetry breaking potential V((f>) = /U 2 </> 2 + A</> 4 . That a condensate can form near the 
horizon without the potential V(cf>) had been previously found analytically for the case of a 
massive scalar field in an asymptotically AdS 4 spacetime [l] . Our numerical solution lends 
strong support to this finding even though in our case the scalar field is massless and the 
spacetime is asymptotically flat. That SSB can occur without the standard potential had also 
been noted in flat spacetime [31. Moreover, as in our case, they obtain that the complex scalar 
field rotates with time. This clearly demonstrates that the phenomena of SSB can occur by the 
properties of the system itself, without an explicit potential, in either flat or curved spacetime. 

The standard stationary BH solutions like Schwarzschild, RN and Kerr do not include a scalar 
field to model uncharged or charged matter. For Schwarzschild, the standard solution yields 
a Ricci tensor and Ricci scalar that is zero everywhere including the interior since this is 
a vacuum solution. However, in spherically symmetric gravitational collapse of uncharged 
matter, the Ricci tensor and Ricci scalar will not be zero inside. In fact, we expect a Ricci 
singularity in the interior where the energy density diverges. The interiors of the standard and 
gravitational collapse solutions clearly do not match; there is an exact match with the standard 
Schwarzschild solution only in the exterior region. In the stationary RN BH case, the standard 
solution does not include a complex scalar field or Dirac field to model the charged matter. It 
is therefore not surprising that the gravitational collapse of a charged scalar field can provide 
something novel in the interior that was not revealed by the standard solution. Again, the 
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numerical solution for gravitational collapse matches the exterior region of the standard RN 
solution exactly. In particular, the Ricci scalar is zero in the exterior region as in the standard 
solution. 

What is of course remarkable is that despite their different interiors, the gravitational collapse 
0013 Mm and standard solutions have the same BH entropy and thermodynamics for 
the same conserved quantities (such as the mass M and charge Q ). As stated in pj3]> this 
is intimately connected to the “no scalar hair” theorem. This theorem however depends on 
having asymptotic flatness. It would therefore be of great interest to investigate numerically 
what happens in the collapse of a charged scalar field to a BH in an asymptotically curved 
background (such as de Sitter or anti-de Sitter space). In particular, one could investigate the 
thermodynamics of such BH’s using the Lagrangian as a probe into the free energy as was 
recently done mm and determine if there is a deviation from standard BH thermodynamics 
due to the presence of the condensate. 
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Figure 1: The radius ro where N = 0 as a function of time. The lapse crosses zero for the first 
time at approximately t=8.4. The value of ro grows with time and approaches a value of 0.344 
at late times. 
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Figure 2: The initial states for the lapse N, the conformal factor i/j and the “electric potential” 
a. Note that the lapse is initially positive everywhere and no black hole has formed yet. 
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Figure 3: Magnitude of charged (complex) scalar field |x(r, t)\ plotted at different times for an 
initial kink shape. Initially, in the interior, the magnitude decreases but then it increases and 
settles to a non-zero constant value at late times. This is in contrast to the uncharged case 
where the (real) scalar field tends to zero in the interior [9]. 
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Figure 4: The Ricci scalar R(r, t ) plotted at different times for charged collapse. It is a non¬ 
zero constant in the interior at late times signaling the presence of a constant energy density 
in the interior. 
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Figure 5: Plot of the real versus imaginary part of the complex scalar field taken in the interior 
at r=0.25 (this is the red curve). Time runs clockwise (from t=0 to t=28 in time intervals 
of 5t = 0.0005). The black curve represents the extrapolated circle constructed from points 
at late times. The red curve deviates from the black curve initially but matches up with the 
black curve at late times (i.e. approaches a circle at late times in the collapse process.) The 
rotation speed u> matches the average value of the electric potential a at late times (see Fig. 0 
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Figure 6: Potential a(r, t) plotted at different times. Inside the horizon, a is a non-zero 
constant. 
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Figure 7: The rotation speed ui of the scalar field and the average value of the electric potential 
a in the interior (r = 0.25) are plotted as a function of time. Note that they approach each 
other with time and match at late times. 
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Figure 8: Testing the WEC. A plot of T tt (r,t ) at different times. Note that T tt is nonnegative 
at all times. 
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Figure 9: A plot of P(r,t ) = T tt /N 2 + T rr /’ip A + 2T rt /(Ni/j 2 ) at different times. See figure 
below for a magnified version. 



Figure 10: Testing the NEC. Note that P is nonengative at all times. 
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